短螺旋型燃烧室旋流结构及气动边界发展规律数值模拟研究

短螺旋型燃烧室的头部轴线沿发动机周向与发动机轴线偏转一定角度,可以在保证燃烧效率不变的条件下有效缩短燃烧室轴向长度。短螺旋型燃烧室流场的最大特点是旋流流动单侧受限。为了研究头部安装角α变化对燃烧室旋流流动特性的影响,基于数值方法对短螺旋型燃烧室进行了计算分析。结果表明:随着α变化,旋流器下游旋涡依次出现对称环状、马蹄状、环状结构;随着α增加,气动边界逐渐出现并抑制旋流的周向扩张,导致流场出现不同的旋涡结构;不同α下的切向角动量随轴向距离增加而衰减,但α为35°和45°时,气动边界在非受限侧出现并对旋流产生约束,角动量衰减变慢;当α为0°,15°,35°,45°时,燃油液滴依次集中分布在旋流器下游两旋涡边缘、侧壁面和头部端壁、非受限旋涡边缘。本文研究了不同α下的旋涡结构及气动边界沿轴向的演化过程,为短螺旋型燃烧室进一步的设计与优化提供基础。     

改进人工蜂群算法及其轴流压气机优化应用

为了探索更优的全局优化算法并将其应用于轴流压气机气动优化上，将标准人工蜂群（ABC）算法中采蜜蜂和观察蜂探索新蜜源的方式进行改进，从而更好地利用总体蜜源的探索信息，得到改进人工蜂群（IABC）算法。经基准函数测试表明，改进算法既能提升全局寻优能力，又能加快收敛速度。采用IABC算法和经过校核的CFD数值方法搭建优化平台，对单级跨声速轴流压气机Stage35进行优化。优化变量为动叶和静叶径向6个截面的弯、掠、前缘重弯值和尾缘重弯值，以流量和压比相对变化保持在0.5%以内为约束条件，以提高绝热效率为优化目标。结果表明：在设计转速下，优化后设计点绝热效率提升0.83%，全工况范围内平均绝热效率提升2.0%，同时喘振裕度提升1.0%，验证了IABC算法在轴流压气机优化中的有效性。     

二元曲面可调进气道流量系数精确预测方法

为了满足二元曲面可调进气道模态转换马赫数范围(来流马赫数为2.2~3.2)的流量要求,针对唇口平移、转动和转动+平移三种调节方案,基于理论分析和基准进气道的流场,提出了一种流量系数精确预测方法,并通过数值计算进行验证且获得了进气道的总体性能。结果表明:调节后的进气道流量系数与预测值完全相等,而且无需多次试算,符合设计预期,可拓展应用于轴对称进气道。相对基准进气道,唇口前移时流量系数和压缩效率同时增加,来流马赫数为2.5时出口总压恢复系数相等而增压比增加了14.6%;在降低相同流量系数条件下,后移唇口使得增压比和压缩效率均降低,来流马赫数为2.5时出口总压恢复系数基本相等而增压比减小了12.9%,转动唇口使增压比进一步减小了9.1%,唇口后移方案性能更优。      

预旋喷嘴径向角度对预旋特性影响的数值研究

为了降低低位进气预旋流路的气动损失，针对带有不同径向角度（0°~30°）预旋喷嘴的预旋系统进行了数值仿真，并对流动特性、温降特性和比熵增特性进行了分析。结果表明:随着预旋喷嘴径向角度的增大，预旋系统无量纲温降先增大后减小，流动阻力减小，预旋系统的流量随之增大。旋转雷诺数为2.3×107时，预旋喷嘴带径向角度的预旋系统无量纲温降比传统喷嘴最大可提高18.3%，存在某一角度使预旋温降特性达到最好。预旋系统内的耗散主要发生在预旋腔和共转盘腔内，径向角度为10°时其比熵增变化量分别占整个预旋系统总体比熵增的42.4%和30.2%；合理设计预旋喷嘴的径向角度，能改善预旋腔内气流的流动效果，并且可以减少整个预旋系统的不可逆损失。      

摆线桨非定常动态失速的三维效应

对摆线桨三维全尺寸模型进行了非定常数值模拟。验证了两种嵌套网格方法的可靠性与适用性,针对二维/三维情况下摆线桨非定常涡流动特性及诱导分离特性开展了数值分析,重点研究了展向气动力随方位角的分布与全局流场的诱导速度分布,并分析了尾迹捕捉精度与涡量耗散特征。结果表明:受摆线桨切向速度、展向诱导速度等三维效应影响,两种计算结果的瞬态气动力极值差值达到44%。三维动态失速涡的产生、脱落与再附明显弱于二维翼型情况,这对于摆线桨非定常气动特性有较大影响。      

扩压叶栅叶尖非定常流动的涡动力学机理

为了揭示压气机叶尖区旋涡结构与流动非定常性之间的关联，采用URANS对一亚声速平面扩压叶栅在不同攻角下的流场进行了求解，并借助Q判据提取了叶尖瞬态涡系结构。结果表明：泄漏涡的破碎现象能够通过诱导新的涡结构间接作用于相邻通道的叶尖流动，是导致叶尖流场失稳的关键因素。在-0.3°和+0.7°攻角下，叶尖泄漏涡发生了螺旋破碎，并伴随有非定常诱导涡的出现，诱导涡对相邻叶片载荷的影响使得叶尖泄漏涡发生周期性摆动；在+1.7°攻角下，泄漏涡破碎会导致反流涡的形成，反流涡的输运会给叶片载荷和来流攻角带来非定常扰动，反过来又会作用于泄漏涡的破碎和反流涡的生成，最终表现为一种自维持的非定常流动现象。     

Numerical investigation on flow nonuniformity-induced hysteresis in scramjet isolator

Numerical simulation and theoretical analysis were conducted to study the hysteresis inside scramjet isolator during the reciprocating process of back pressure variation. It is revealed that only a regular reflection is theoretically possible for two leading shocks when the inflow Mach number is greater than 2.0, and no hysteresis can occur in the transition between shock reflection types. Nevertheless, wall suction, gas injection, and background waves cause non-uniformity of the incoming flow and would make hysteresis possible. Besides the classical hysteresis in the transition between shock reflection, new kinds of hysteresis were found in both the deflection angle of separated boundary layer and the location of the shock train. Moreover, the occurrence of hysteresis in the deflection angle of the separated boundary layer is accompanied with the shock reflection hysteresis. In the case with background waves or gas injection, hysteresis in the starting position of leading shock was observed too. As back pressure decreases, the leading shock does not follow the same path as that as the back pressure increases, and it is anchored at the location where the background shock or the injection interacts with the leading shock. It is inferred that, if two strong adverse pressure gradient regions move towards and interact with each other, hysteresis will take place when they start to separate.     

Unsteady wakes-secondary flow interactions in a high-lift low-pressure turbine cascade

Detailed experimental measurements were conducted to study the interactions between incoming wakes and endwall secondary flow in a high-lift Low-Pressure Turbine (LPT) cascade. All of the measurements were conducted in both the presence and absence of incoming wakes, and numerical analysis was performed to elucidate the flow mechanism. With increasing Reynolds number, the influence of the incoming wakes on suppressing the secondary flow gradually increased owing to the greater influence of incoming wakes on reducing the negative incidence angle at higher Reynolds numbers, leading to a lower blade loading near the leading edge and suppression of the Pressure Side (PS) leg of the horseshoe vortex. However, the effect of unsteady wakes on suppressing the profile losses gradually became weaker owing to the reduced size of the Suction Side (SS) separation bubble and increased mixing loss in the free-flow region at high Reynolds numbers. Incoming wakes clearly improved the aerodynamic performance of the low-pressure turbine cascade at low Reynolds numbers of 25,000 and 50,000. In contrast, at the high Reynolds number of 100,000, the profile loss at the midspan and mass-averaged total losses downstream of the cascade were higher in the presence of wakes than in the absence of wakes, and the unsteady wakes exerted a negative influence on the aerodynamic performance of the LPT cascade.     

光敏推进剂激光烧蚀过程建模与仿真

相较于传统大卫星,微小卫星具有结构紧凑、质量轻便和成本低廉的特点。然而,受功率和质量负载的限制,微小卫星一般不装备推进系统,其航线也局限于近地轨道。为扩展微小卫星的功能,满足日益复杂的任务需求,需给其配备合适的微推进系统。固体推进系统具有结构简单、寿命长、可靠性高的优点,但无法重复启动。为得到可重复启动的固体微推进系统,设计了一种非自持燃烧的光敏推进剂,采用激光控制其燃烧。在背压为大气压的环境下,利用高速摄像机拍摄燃烧过程并记录燃速。之后,对光敏推进剂的激光烧蚀过程进行建模。分析结果表明:激光可控制光敏推进剂的燃烧,燃速与激光强度成线性关系;该光敏推进剂的最小激光点火强度为0.28 W/mm~2;燃速计算值与实测值的误差在10%以内,证明该数学模型具备工程应用价值。     

场协同原理在高超声速化学非平衡流动中的推广

应用理论分析与数值模拟方法，将对流传热的场协同原理从不可压缩流动推广至高超声速化学非平衡流动中。结果表明，高超声速化学非平衡层流与湍流的热流密度取决于流动的当地单位体积的动量与单位质量总焓梯度的协同。用当地单位体积的动量与单位质量总焓梯度的协同研究高超声速化学非平衡流动的壁面传热问题，对层流流动下的对流传热，不但计及了高超声速化学非平衡流的密度变化对热流密度的影响，而且包括了静焓梯度、压力梯度、边界层内的分子黏性剪切效应对热流密度的作用；对湍流问题，除了上述层流流动各项对热流密度的影响外，还计及了雷诺剪切应力对热流密度的作用。考虑到高超声速化学非平衡流静焓的定义，高超声速化学非平衡层流及湍流的场协同同时计及所有组分的平动能、转动能、振动能及电子能等梯度的贡献。